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Abstract 

The recent observation of a Higgs boson at the LHC and experimental confirmation of the nonvanishing 
neutrino-mixing parameter sin^is offer important means to test physics beyond the standard model. We 
explore this within the context of the scotogenic model, in which neutrinos acquire mass radiatively via 
one- loop interactions with dark matter. Starting with a two-parameter neutrino- mixing matrix which is 
consistent with the latest neutrino-oscillation data at the one-sigma level, we derive different sets of solutions 
for the Yukawa couplings of the nonstandard particles in the model and use the results to consider the 
Higgs decays into final states involving the new particles. Assuming that the lightest one of them serves as 
fermionic cold dark matter, we show that such decays are allowed by various experimental and theoretical 
constraints to have substantial rates that are already restricted by the current LHC data. We also look 
at their correlations with the Higgs decays into 77 and ^yZ. Upcoming LHC measurements of the Higgs 
boson can therefore either detect scotogenic signals or place further constraints on the model. 
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I. INTRODUCTION 



The recent discovery of a Higgs boson with mass in the 125-126 GeV range at the LHC [1] 
and experimental confirmation of the neutrino-mixing parameter sin ^13 that is nonnegligible [2] 
undoubtedly have far-reaching implications for efforts to identify the nature of physics beyond 
the standard model (SM). Any realistic scenario for new physics would need to incorporate such 
a particle and account for neutrino masses and mixing angles of the right amount. In addition, since 
about 22% of the total cosmic energy density has been inferred from astronomical observations to 
be attributable to matter that is nonluminous and nonabsorbing [3], the desired new- physics model 
should also possess at least one candidate for dark matter (DM). 

Among the simplest possibilities accommodating the necessary ingredients is the scotogenic 
model proposed by Ma [4] , in which neutrinos get mass radiatively via their one-loop interactions 
with new particles comprising scalars and fermions, at least one of which plays the role of DM. Here 
we explore some implications of the aforementioned experimental findings within the context of the 
minimal version of this model. In particular, identifying the newly observed Higgs boson with the 
Higgs boson in the model, hereafter denoted by h, and assuming the lightest one of the new fermions 
to be the DM, we consider the decays of h into final states containing the nonstandard particles 
in the model. Moreover, since their Yukawa couplings depend on the neutrino-mixing angles and 
contribute to the decay amplitudes, it is important to adopt parametrization for the couplings that 
takes into account the fact that sin ^13 is not negligibly small. Subsequently, we will demonstrate 
that such exotic decays of h axe allowed by various experimental and theoretical constraints to have 
substantial rates that are already bounded by the existing LHC data. As upcoming measurements 
at the LHC will pin down the various properties of h with increasing precision, the acquired data 
will then either reveal hints of the new particles or probe the model more stringently. 

The organization of the paper is as follows. The next section gives a description of the relevant 
Lagrangians and the neutrino mass formula. In Section III, we write down a neutrino-mixing 
matrix that depends on only two parameters and is consistent with the neutrino-oscillation data, 
including the measured nonzero ^^13, at the one-sigma level. From the resulting neutrino- mass 
matrix, we derive solutions for the Yukawa couplings of the new particles in the model. In Section IV, 
we look at a number of experimental and theoretical constraints on their couplings and masses. 
Specifically, there are low-energy measurements that can limit the Yukawa couplings. Furthermore, 
the Yukawa couphngs belonging to the DM candidate also have to be compatible with the observed 
relic abundance. In Section V, with the parameter values satisfying the preceding constraints, we 
investigate the Higgs decays into final states involving the nonstandard particles and take into 
account extra constraints from the latest LHC data, including those on the Higgs decays into 77 
and 7Z. We also examine the scotogenic impact on the correlations between these different decays. 
We conclude with a summary of our results in Section VI. 

II. INTERACTIONS 

In its simplest version, the scotogenic model extends the minimal SM with the addition of only 
a scalar doublet, 77, and three neutral singlet fermions, N^, all of which are odd under an exactly 
conserved Z2 symmetry [4]. The SM particles are all even under this symmetry. Accordingly, the 
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lightest one of the new particles is stable and can act as DM. In this study, we consider the case 
that Ni is a good candidate for cold DM [5].^ 

The Lagrangian for the interactions of the scalar particles in this model with each other and the 
standard SU(2)LxU(l)y gauge bosons, Wp and Bp, has the form 

jC = {VP^yVp^ + {VPriyVpT] - V , (1) 

where T>p^dp + {i/2)gT-Wp + igyQyBp, 

+ A3($t$)(^t^) + A4($t^)(77t$) + iA5[($t^)2 + (77t$)2] , 
and, after electroweak symmetry breaking. 



(2) 



H+ 



(3) 



with g (gy) being the usual SU(2)l (U(l)y) gauge coupling constant, r = {T-^^,T2,r.^) the Pauli 
matrices, Qy the hypercharge operator, and v the vacuum expectation value (VEV) of $. The 
VEV of 77 is zero due to the Z2 symmetry. The masses of 5, V, and H"^ are then, respectively. 



m< 



fil + |(A3 + A4 + As)!;' 



m: 



V 



m 



H 



/i2 + lAsW^ . 



(4) 



In our numerical analysis in Sections IV and V, we will make the usual assumption [5] that A5 is 

lAsI^^ <C m| ~ m^. From Eq. (1), 



very small, IA5I ^ IA3 + A4I, which implies that |m| 
the couplings of rj to h, the photon A, and the Z boson are described by 

h 

.2 , -^/(i-24) 



+ ie {H+ dPR- - H' dPH+)Ap + H+H'A^ 



H+H-A^Z^ 



2c, 



(5) 



where only terms relevant to our processes of interest are on display, e = gs^ > is the electro- 
magnetic charge, and = 



'1 — = cos^w with the Weinberg angle 
The new singlet fermions A^^^ are permitted to have Majorana masses and interact with other 
particles according to 



j.u^(S-tV)\Pj,N, + H.c. , (6) 

where j, k — 1,2,3 are summed over, the superscript c refers to charge conjugation, = 1(1 + 75), 
and -^1,2,3 = e, /i, r. Hence, writing the Yukawa couplings 3^^^. — Y^.k, we have 

^ Yel Yf,2 Yf,^ ^ 
Yp,! Yp2 Yps 

\ Yrl Yt-2 Yt-3 J 



y = 



(7) 



The possibility of Nk being warm DM has also been proposed in the literature [6]. 
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In this model the hght neutrinos get mass radiatively through one-loop diagrams involving in- 
ternal S or V and N^. The mass eigenvalues rrij of the neutrinos are given by [4] 

diag(mi, mg, mg) = U^M^U , (8) 
= 3^diag(Ai,A2,A3)3^^ , 



A. 



Xix) 



X 



lux 



1 — X il — x) 



(9) 



2ml = ml + rrip , (10) 



where U is the PMNS (Pontecorvo-Maki-Nakagawa-Sakata [7]) unitary matrix and the expression 



for Aj^. is valid for ttIq ~ ~ m-p. 



III. MIXING AND YUKAWA MATRICES 



We express the PMNS mixing matrix as a product of two matrices involving only two mixing 
angles, 9 and q, respectively, with the latter matrix also containing a CP-violating phase S. Thus 

/ cos6' sin6' ^ \ f cosq e^'^sin^\ 

1 



U 



y sin ^ ^ cos 9 



^ ^ 73 / V "^""^ ^™ ^ ^ cos ^ / 

/ "\/2 cos 9 cos q \/2 sin 9 \/2 e"^ cos ^ sin ^ \ 

— sin ^ cos — e"*'' sin cos^ — e*"^ sin^ sin^ -|- cos^- 

\ sin 9 cos — e~"^ sin — cos 9 e"^ sin 9 sin <j -|- cos <j / 

The form of U with sin 9 = 1/ \/3 was discussed in Rcf. [8] , whereas the <^ = case was treated 
in Ref . [9] . Both of these possibilities for U arc no longer compatible with the most recent findings, 
especially that sin ^13 is not neghgibly small [2] . Therefore, we will instead take 

cos^^ sin<; = sin^^-|^3 , 9 9^2 1 (12) 

which lead numerically to elements of U consistent with their empirical counterparts within one 
sigma. For simplicity, we also fix e*"^ = 1 in accordance with the value 5 = (SOOl^gg)" from the 
latest fit to the global data [10]. 

Incorporating Eq. (11) into the matrix diagonalization relation in Eq. (8), we then derive the 
mass eigenvalues 

+ \ [{Y^k - Yrk)sec^ + {Y^k + i;fc)sj'}Afc , 



1 

7^ 



(11) 



mi 



rrir, 



A, 



Ca 



(13) 
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where we have imphcitly summed over k — 1, 2, 3 and adopted the notation 

Cq — cos 9 , Sg = sm9 , = cos^ , = sin^ . 



(14) 



As will be seen shortly, these mass formulas can be rendered much simpler using the relations among 
which have to fulfill the required vanishing of the off-diagonal matrix elements on the right-hand 
side of Eq. (8). Thus we arrive at the diagonalization conditions 







+ {Y^, - Y^\)co c^ + V2 Fefe {Y^k + Yrk)so cj Afe , 



+ 



- Y^, - Y,\)4 - V8 {Y,, - Y^,)c,s, - 2Y^,Y,, (l + s'^) 
V2 Y,k - {Y^k - Yrk)sJ^ {Y^k + Yrk) {4 - 4) , 



(15) 



summation over k — 1,2,3 being again implied. It turns out that these equations are exactly 
solvable for Yek and 1^^ in terms of Yk = l^fe. As sketched in Appendix A, there are twenty-seven 
possible sets of solutions to Eq. (15), but three of the sets can each produce only one nonzero mass 
out of m^2 3 Eq. (13), whereas another eighteen (six) of the sets can each lead to two (three) 
nonzero masses. 

It is worth pointing out that the form of Eq. (9) also appears in some other models of radiative 
neutrino mass, which may be generated by one-loop [11] or two-loop [12] diagrams. Hence these 
solutions for Yg^ plus the resulting masses mi,2,3 are also applicable to such models, with Ai^2,3 
encoding the model specifics. 

The solutions in one of the eighteen sets that can each yield two nonzero masses are 



Y = 



Y 



V2 c, Yi 



Sg + 

— Sq 

Sg + 



i = 1,2 



Ye3 = 



-V2sgYs 



Yi, 



Y^s — 



-Y. . 



(16) 



These lead to the masses 
nil — , 



rrio — 



2A3 Yi 



rrin 



2{A,Y,' + A,Yi) 
{sgs^ + c^ 



(17) 



Setting <j = in the last two equations, one recovers the corresponding expressions derived 
in Ref. [9]. Since it is now known experimentally that sin^^i3 — 0.0227lo;oo24 [10] > which is not 
very small, the = limit is no longer a good approximation. Particularly, as shown below, this 
nonzero ^13 corresponds to Fgi — 0.24 and Y^j ~ 0.82 l^j in Eq. (16), compared to Y^i — and 
Yfj,i — Y^i in the ^ = case [9]. 
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Hereafter, we will employ nonzero q according to Eq. (12) and focus on one of the solution sets 
as an example that gives three nonzero masses. The solutions in this set are 

_ V2c,c^Y^ _ -V2sgY2 _ V2c,s^Ys 

-'el — — , J^e2 — , J^e3 — — , 

~ + 

1^1 = ^i±^y, , Y,, = -Y, , y,3 = ^^^^^3 , (18) 



which yield 



2AiFi2 2A2Yi 
m, = , m, = , m, = . (19) 

These expressions for 2,3 would permit cancellations among the terms in Am^j = — mf with 
larger Y^ than would the masses in Eq. (17). 

Numerically, we adopt for definiteness 



cos 



^sin^ = VO.0227 , = 32.89°, (20) 



which translate into elements of U that are well within the one-sigma ranges of their experimental 
values. These choices imply the numbers collected in Table I for the Yukawa couplings and neutrino 
masses given in the previous two paragraphs in terms of Yj^ and A^F^^. After applying Eqs. (18) 
and (19) to the constraints to be discussed in the following section, we will employ the allowed 
parameter values to explore several decay modes of the Higgs boson h. 



TABLE I: Numerical values of Y^^ = Y^j^/Yf. and neutrino masses in terms of = Aj^Y^ from the 
formulas in Eqs. (16)-(19), with the input parameters in Eq. (20). 



Equations 


Yel 


Y^.1 


Ye2 Y^2 


Ye3 


Yf^3 


nil 


m2 




(16), (17) 


0.197 


0.820 


0.197 0.820 


-0.915 


-1 





2.84 A3 


1.7l(Ai +A2) 


(18), (19) 


3.293 


-2.011 


-0.915 -1 


0.197 


0.820 


15.89 Ai 


2.84 A2 


1.71 A3 



IV. CONSTRAINTS 



The couplings and masses of the nonstandard particles in the scotogenic model are subject to 
various constraints. Theoretically, there are a number of restrictions on the parameters Aj in the 
potential in Eq. (2) and on the Yukawa couplings yjk- Vacuum stability is ensured by demanding 
that Ai^2 > 0, A3 > — v^AiA2, and A3 + A4 ± IA5I > — -\/AiA2 [13]. The condition of perturbativity 
translates into |Aj| < 87r and \yjk\ < V^vr [14]. There are additional requirements from unitarity 
on some combinations or functions of Aj [15, 16]. 

Experimentally, there are constraints on the masses of the new scalars. The data on W and Z 
widths and the null results of direct searches for new particles at e+e~ colliders imply that [16-18] 

rriH + mg -p > , rrijj > 70 GeV , rrig + m-p > niz , (21) 
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where the last inequahty is vahd for jm^ — m^l < 8GeV [18], which pertains to our assumption 
that ~ m-p. Accordingly, in our numerical analysis later on we will consider the mass regions 
50 GeV < 

'^sv — 120 GeV and 70 GeV < — 120 GeV. These choices respect the limits on the 
oblique parameters S and T at 90% CL (Gonfidence Level) [3, 19]. 

There are also experimental constraints on the Yukawa couplings y^^^ from several low-energy 
measurements and the observed DM relic abundance, which we address in the rest of this section. 
The latest LHC data on the Higgs boson imply extra restrictions on the model, which we will take 
into account in Section V. 



A. Low-energy observables 



Neutrino oscillation measurements determine the differences A^,- 



fit to the data [10] 



2 2 



Prom the latest 



^21,exp 



(7.50lO:J8) X 10-5 eV^ , 



A^ 

jl,exp 



(2.4731°:°™) X 10"' eV^ 



(22) 



We have here assumed the normal ordering of the neutrino masses, which is preferred to the inverted 
ordering by the solutions in Eq. (18). In restricting the new couplings, we will then impose 



31.0 < 



A 



31 



A2 

^21 



< 35.0 



(23) 



based on the 90% CL ranges of the numbers in Eq. (22). 

For the masses, there are also constraints on the effective mass parameters (^^) and (^^^) 
from beta decay and neutrinoless double-beta decay experiments, respectively, and on the sum of 
masses, Sj^.mj., from astrophysical and cosmological observations. The hmits are [20] 



<m^> = ^Jl2\Uik\^rnl < 2.1 eV , (m^^> 



< 0.25 eV , 



(24) 



and Yij^m^, < (0.5-1.5) eV. As it will turn out, these are less restrictive on the Yukawa couplings 
than the other constraints described in this subsection. 

The interactions of and with charged leptons give rise to the flavor- changing radiative 
decay £j — > 7 at one- loop order. Such decays have been searched for, with negative results so far, 
including the fresh one for — > 67 reported by the MEG Collaboration [21]. The experimental 
bounds on their branching ratios are [3, 21] 



B{n e7)exp < 5.7 x 10 



-13 



B{r ^ e7)exp < 3.3 x 10"' , 



B{t ^ /X7)exp < 4.4 X 10" 
at 90% CL. Hence they put a cap on the prediction [5, 22] 



Y.yiky*jk^{Mi/my) 



(25) 



(26) 
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where Gp — v ^/-\/2 is the Fermi constant, 

" = 4^ ' -^^^^ = eo^^p ' ('^^ 

and numerically for B{£j diPu) we will use the central values of their data: B{ii ei'u)exp — 1, 
B{r ei/z/)e^p = 0.1783 ± 0.0004, and B(t /xz/z>)exp = 0.1741 ± 0.0004 [3]. 

At one-loop level, the presence of H"^ and Nj also induces a modification to the anomalous 
magnetic moment a^. of lepton given by [22] 

The existing data on and the charged-lcpton masses imply that only can be significantly 

restrictive on the potential scotogenic effects at present. Its most up-to-date SM and experimental 
values differ by nearly three sigmas, a^^'P — = (249 ± 87) x 10~^^ [23]. Accordingly, in view of 
the negative sign in Eq. (28), we may require 

|Aa^| < 9 X 10-^° . (29) 

This will turn out to be complementary to Eqs. (23) and (25) in restraining the Yukawa parameters. 



B. Fermionic dark matter 



Since we have picked A^i to be the lightest of the nonstandard particles and serve as cold DM, it 
needs to account for the observed cosmic relic abundance, which therefore imposes bounds on J'^^. 
The Ni annihilation cross-section a^^^ is related to its relic density fl hj [5, 24] 



1.07 X 10^ Xf GeV"^ 



In 



0.0955 {a + 6b/xf)M^mp^ 



(30) 



where h denotes the Hubble parameter, mpj = 1.22 x 10^^ GeV is the Planck mass, is the number 
of relativistic degrees of freedom below the freeze-out temperature = M-^/xj:, and a and b are 
defined by the expansion CTg^n'^rei = « + ^'^^rei + ^('"rei) terms of the relative speed of the NiNi 
pair in their center-of-mass frame. 

The leading contributions to cx^^j^ are the tree-level processes N-^N-^ — > and N^N^ I'^i'j 

via exchanges of H"^ and [S,V), respectively, each of which proceeds from diagrams in the t and u 
channels because of the Majorana nature of the external neutral fermions. We collect the expressions 
for their squared amplitudes in Appendix B for ttIq c::^ m-p. If all the final-lepton masses are 

negligible, their combined cross-section times v^^y is 



E 

ij=l,2,3 



487r 



Mf 



m 



H 



Mf + m: 



(31) 



to second order in which implies a = and b 0. The = ttIq limit of this formula agrees 
with that found in Ref. [5]. In our numerical computation with more general masses, we employ the 
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cross section obtained from the squared amplitudes in Appendix B and do not neglect the charged 
lepton masses, in which case a is also nonvanishing. 

Since the solutions for derived in Section III are all proportional to Yi, it is the only coupling 
relevant to the relic density of Ni. To extract \Yi\ from the empirical value of ^2, one can utilize 
Eq. (30) once the mass parameters ttIq ^ and are specified. Thus, we present in Figure 1(a) some 
samples of the values of iFil consistent with the 90% CL range of the data Vth? = 0.111 ± 0.006 [3] 
over 5 GeV < Mi < 50 GeV for the solutions in Eq. (18) and various sets of Tn^^. To offer 
a different perspective on the allowed values, in Figure 1(b) we display \Yi\ versus uiq for several 
choices of Mi and m^. The plot of |Yi| versus could be roughly inferred from Figure 1(b) by 
interchanging uIq and m^, especially for -C Mi. Evidently, the demand that A^^i be the leading 
candidate for cold DM over the mass regions of interest can always be met by some values of Yi, 
but their ranges are fairly limited. 

Direct-search experiments for DM, which look for signals of it colliding with nuclei, may lead 
to further restraints on Yi, but the existing bounds are still too weak. Since A^^i can scatter off 
a nucleon mainly via its one-loop Z-mediated axial- vector interactions with quarks [9, 25], the 
process is characterized by a spin-dependent cross-section that is relatively suppressed and in our 
case does not reach 10"''^ cm^. This is more than an order of magnitude below the strictest limit to 
date, measured by the XENONIOO Collaboration [26]. 




20 30 
Ml (GeV) 



80 90 
OTq (GeV) 



FIG. 1: Magnitude of Yi belonging to Yukawa couplings in Eq. (18) versus (a) Mi and (b) ttiq satisfying the 
relic density constraint 0.101 < Qh? < 0.121 for some choices of {mQ,m^^ and (M^,mj:^), respectively. 



V. IMPLICATIONS FOR HIGGS BOSON DECAY 



As experimental work on the Higgs boson proceeds at the LHC with increasing precision, the 
accumulated data will reveal how much the properties of the particle may deviate from SM expec- 
tations. The information gained will then serve to test in particular various scenarios in which new 
physics can induce nonstandard decay modes of the Higgs and/or significant modifications to its 
SM decay channels [27]. In the scotogenic model, such effects can arise at tree and loop levels, to 
which we now turn. Using the parameter space allowed by the constraints discussed above, we first 
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look at Higgs decays into final states containing the new particles. As mentioned earlier, in our 
numerical work below we assume 50GeV < ttiq < 120 GeV and 70 GeV < < 120 GeV. 

Since collider data imply 2mQ > m^, our mass ranges of interest include < 2mQ < m^, in 
which case the decay channels h — > SS,VV are open and may be important [17]. From Eq. (5), 
we obtain the amplitudes for these modes at tree level to be 

2 (ml — u,T\ , , 

M^^ss ^ My^^vv ^ ^ \ ^ ■ (32) 

Their combined rate is 



2 



Y{h ^ r?V) = r(/i ^ SS) + V{h ^ VV) ^ /i _ . (33) 

The decay products would be invisible for S and V lighter than Nk- Otherwise, S and V will 
decay into i^jN)^ and, if kinematically possible, will subsequently decay into t^l'J,Ny {p^u-iNy^ 
through (5, V) exchange. 

In the mass range m^/2 < tUq < nif^, these two-body decays no longer happen, and so the 
dominant modes with S and V in the final states are the three-body decays h {S,V)iyjNf, if 
kinematically permitted. Prom Eqs. (5) and (6), we obtain the tree- level amplitude 

Mn-,s(r).,N, = {-^)- , f ' \^ ^^^^ 

and an analogous expression for h — > S{V)PjNf.. They contribute to the combined rate 

T{h^v"i^N) = J2 E [^{h^fli^jN,)+r{h^fju^N,)] . (35) 

fi=s,r i,fc=i,2,3 

For decays with in the final states, since collider data restrain their mass to be m^j > 70 GeV, 
we only have the three-body modes h — )■ H~^i~Nj, and their charge-conjugated counterparts, pro- 
vided that mjj + m^ + m^ < m^^. As in the neutral case, we find 



2 iml-i,l)yjJ^P^ N, 
and similarly for h H't^N^,. They lead to the rate 



V{h ^ Hm) = [T{h^ H^i-Nk) +r{h^ H-i^N,)] , (37) 

j,fe=l,2,3 

which is identical in form to Eq. (35). 

After fixing Yi to the range fulfilling the relic density requirement as in Section IV B, we scan all 
the relevant parameter space of the model subject to the constraints described in Section IV A. We 
find that, for the allowed couphng values, these two- and three-body decay channels of the Higgs can 
have enlarged rates. We illustrate this in Table II for different sets of m^ jj, //2, and Mi^2,3, employing 
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TABLE II: Sample values of mass parameters rrigfj, 112, and Mi 2,3, all in GeV, and Yukawa constants 
^1,2,3 satisfying the constraints discussed in Section IV. The last three columns contain the resulting 
branching ratio Bsvh, in percent, of Higgs decay into final states containing S,V, or and ratio TZ^yo 
of r{h 7V°) to its SM value for V" = 7, Z. 



mo 






Ml 


M2 


M3 




Y2 


^3 


^SVH 






50 


70 


46 (47) 


9 


13 


63 


0.155 


0.372 


0.632 


21 (14) 


0.89 (0.89) 


0.95 (0.95) 


60 


80 


54 (56) 


10 


15 


71 


0.172 


0.413 


0.700 


26 (14) 


0.91 (0.92) 


0.96 (0.97) 


70 


70 


40i (25) 


14 


20 


80 


0.155 


0.378 


0.677 


23 (11) 


0.75 (0.83) 


0.88 (0.92) 


70 


80 


111 (99) 


12 


17 


78 


0.180 


0.437 


0.730 


24 (12) 


1.15 (1.09) 


1.06 (1.04) 


80 


80 


159 (88i) 


15 


22 


82 


0.175 


0.422 


0.760 


22 (13) 


1.53 (0.68) 


1.21 (0.86) 


120 


70 


123 (111) 


20 


29 


85 


0.157 


0.386 


0.715 


20 (12) 


1.48 (1.34) 


1.20 (1.14) 



the Higgs mass = 125.5 GeV, compatible with the most recent measurements [28], and SM Higgs 
total width — 4.14 MeV [29]. The numbers wc have selected for uiq Mi, and Yi can be seen 
to correspond to some of the points inside the colored bands on one or both graphs in Figure 1. In 
the tenth column of this table, we list the branching ratio Bsvh = ^svh /{Tf^ + TsvH) involving the 
rate Tsvh which is the sum of V{h — >■ rf'r]") or T{h rfvN), depending on mg, and V{h — )■ HiN). 
The two numbers on each line in the Bsvh column correspond to the two numbers on the same line 
in the //g column, where we have included the possibility that //| can be negative [13]. Evidently, 
Bsvh can be readily altered by only varying fj,2, with the other parameters fixed. 

The substantial values of Bsvh in this table are made possible by partial cancellations in Am^j 
between the terms proportional to Y^^^ according to Eq. (19) and in the ij — )> £^7 rates between the 
different Yi-kY^.i; terms with opposite signs, as Table I indicates, which allow ^1,2,3 not to be too 
suppressed by the stringent constraints. Some of such cancellations would not happen with solution 
sets of the type in Eq. (16). We should mention, however, that the rather sizable 1^1,2,3 can result 
only with some degree of fine-tuning, roughly at the per-mill level. For such Y^ and the masses in 
Table II, we can reproduce the measured A21 31 with A5 ~ 2 x 10~^^. 

To make comparison with LHC results, we find that the larger (unbracketed) numbers for Bsvh 
in Table II have begun to be probed by bounds inferred from the latest Higgs measurements. 
According to several analyses [30], the current data imply that the branching ratio of nonstandard 
decays of the Higgs into invisible or undetected final-states can be as high as 22% at 95% CL if 
the Higgs production mechanism is SM-like, which is the situation in the scotogenic model. These 
limits are not yet very strict and, as the bracketed Bsvh numbers indicate, can be easily evaded 
by changing /i2, which still has a wide range of viability. Therefore, the availability of other decay 
modes which may provide complementary constraints would be highly desirable. 

The impact of the new particles in the model can also be generated through loop diagrams.^ Of 
great interest are their contributions to standard decay channels of the Higgs that are already under 



^ The new particles can contribute to .Z-boson decays into neutrinos (charged leptons) via one-loop diagrams with 
internal S or V (H^) and N/.- At tree level, the Z boson can also decay into three-body final states containing 
them, similarly to the three-body Higgs decays above. We have checked that these scotogenic effects on the Z 
decay are not significant. 
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investigation at the LHC. Here we look at /i ^ 77 and h 7Z which arise in the SM mainly from 
top-quark- and W^-boson-loop diagrams and also receive one-loop contributions from H^. These 
transitions are the same as those in the inert doublet model [16, 17, 31]. Based on general results 
in the literature [32, 33], the predicted rates of these modes are 



r(/i^77) = 



128^2 7r3 



m 



H 



(38) 



64:%'^ m\ 



3c 



(39) 



where the expressions for the form factors ^'^^ available from Ref. [33], the A^^''^^ terms 

originate exclusively from the contributions, Hx = 4m^/m^, and Xx — 4m\/m\. It is worth 
noting that in the = limit the amplitude for h — > 7Z would reduce to that for /i — > 77 
modulo the different 7 and Z couphngs to fermions, W bosons, and H^. 

In Table II we have also listed the resulting numbers for the ratio 

where r(/i — >■ 7V°)gM is the SM rate, without the part. The examples in this table demonstrate 
that the scotogenic effects on r(/i — ?> 77) and r(/i — )> 7Z) have a positive correlation. Comparing 
the BsvH and TZ-y-y^-yz numbers, we see that the latter have the milder dependence on /i2, unless 
changes signs. Moreover, there does not appear to be a clear correlation between Bsj,^ and the 
impact of on r(h — > 7V°), which is partly due to the fact that the Yukawa parameters which 
are present in T^vh do not contribute to r{h — >■ ^V°). Since /i — )■ 77 has been detected, unlike the 
'-fZ mode [34] , we can already compare our examples with the data. The latest measurement of the 
signal strength for /i — )> 77 performed by the ATLAS Collaboration is cr/a^^ = 1.6 ± 0.3 [35]. On 
the other hand, for the same mode the CMS Collaboration has found cr/ugj^ = 0.78 ± 0.27 and 
1.11±0.31 using two different methods [36]. While awaiting an experimental consensus on this decay 
channel, we can say that all of the TZj^ numbers in Table II are still compatible with one or more of 
these LHC results, but the situation will likely change when more data become available. Wc may 
also expect that supplementary information will be supplied by future observations of /i — )■ jZ. 

We have seen that there are some decay modes of the Higgs boson that can be employed to 
test the scotogenic model in complementary ways. At this point the restrictions inferred from the 
existing data on the Higgs decays are not yet strong, but already start to probe the parameter space 
allowed by other data. In addition, with the fairly large Yukawa couplings which we have obtained 
and the relatively light charged scalars, direct searches at the LHC may offer extra tests [37] on the 
scenario treated here. 
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VI. CONCLUSIONS 



We have explored for the scotogenic model of radiative neutrino mass some implications of the 
recent discovery of a Higgs boson at the LHC and experimental determination of sin ^13 that is 
not very small. Employing a two-parameter neutrino-mixing matrix which is consistent with the 
latest neutrino-oscillation data within one sigma, we derive solutions for the Yukawa couplings of 
the nonstandard particles in the model, which consist of scalars and fermions. Such solutions are 
also applicable to some other models of radiative neutrino mass. We select one of the new fermions 
to be the lightest of the nonstandard particles which plays the role of cold DM. Subsequently, taking 
into account various constraints, including those from low-energy measurements and the observed 
rehc density, we use the solutions for the Yukawa couphngs to consider Higgs decays into final states 
containing the new particles. We find that within the allowed parameter regions the rates of such 
decays can be significant, which can already be probed with the latest Higgs measurements. We 
also examine how these exotic decay channels may correlate with the scotogenic effects on the Higgs 
decays into 77 and 7Z, which are under intensive study at the LHC. Consequently, upcoming Higgs 
data with improved precision from the LHC, or a future Higgs factory, can be expected to reveal 
hints of the new particles or impose further restrictions on the model. 
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Appendix A: Solutions for Yukawa couplings 

One can solve the diagonalization conditions in Eq. (15) exactly for the three pairs of Yukawa 
couplings {Yek, l^fc), k — 1,2, 3, in terms of = Yrk- There is more than one set of the solutions. 
In each set, we can express the pairs as (Fgfc) ^^fe) ~ (^z' J^z)^ki where z = a,h, or c and 



S0 Cq S0 + 

K - '-^^^ , -I, • (Ai) 

Sq + 

Since two or all three of the pairs may share the same z, such as in Eq. (16), there are altogether 27 
sets of the solutions to Eq. (15). Not all of them can lead to at least two nonzero masses among the 
eigenvalues in Eq. (13). Three of the sets can each only yield one nonzero mass, whereas 18 (six) of 
the others can lead to two (three) nonzero masses. 



Appendix B: Dark matter annihilation amplitudes 

The Majorana nature of Nj imphes that the process NkNi — >■ arises from if-mediated t- 

and li-channel diagrams. We find the absolute square of its amplitude, averaged (summed) over 
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initial (final) spins, to be 

{M^+mj^-t){Mf + mj-t) 



•^NkNi-^eref\ - \yikyji\ 



+ \yiiyjk\'' 



{Ml + m\ - u) {Mf + mj. - u) 



where 



+ My:ky,iyuy;k) ;w ^ v , (bi) 



« = {PNk+PNiY , t = (pjv, -P^i)^, u = {pNk-PijY ■ (B2) 



In the case of N^Ni — > i/jp^, proceeding from (7^, <S)-mediated t- and w-channel diagrams, we 
need to take into account the Majorana nature of the final neutrinos as well. It follows that for 
uIq ~ ~ m-p and negligible u masses 



•^ATfcATi-^i/ipj — lyikyjii"^ - — —. — p — 72 — - + \yiiyjk 



{Mi-t){M^-t) , „(M|-«)(Mf-«) 



4(mo-t) A{ml-u) 

Re{y:,y,,yuy;k) — v , (b3) 



where now t = [pNk - Pv,)^ and u = [pn^ ~ P^jY ■ 
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